TASEP on a ring with internal degrees of freedom 
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A totally asymmetric exclusion process on a ring with v non-conserved internal degrees of freedom, 
where particles hop forward with a rate that depends on their internal state, has been studied. We 
show, using a mapping of the model to a zero range process with v different kinds of boxes, that 
steady state weights can be written in a matrix product form and calculate the spatial correlations 
exactly. A comparison of the model with an equivalent conserved system reveals that unequal 
hopping rates of particles belonging to different internal states is responsible for the non-trivial 
correlations. 
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I. INTRODUCTION 



Driven diffusive systems have been studied exten- 
sively in recent past for their unorthodox non-equilibrium 
properties [1|. Totally asymmetric simple exclusion pro- 
cess (TASEP) is one such system, initially introduced @ 
as a basic model for transport [3j, which has found ap- 
plications in wide areas of physics and biology The 
steady state weights of TASEP, which has been calcu- 
lated exactly 0-0] , show collective behaviour like bound- 
ary driven phase transitions, non-trivial density and 
shock-profiles. Several variations of this exclusion model 
including non-conservation of particles, multiple species 
[1, Q have been studied in different contexts. 

Recently TASEP with non-conservation of parti- 
cles having internal degrees of freedom has also been 
studied[10]. In these models particles are allowed to enter 
or exit the system from both boundaries and their hop 
rate in the bulk depend on the internal degree they pos- 
sess. The explicit conditions on the rates for the steady 
state of a parity-time invariant system to have a fac- 
torized form has been derived. Analytical results for the 
steady state of these models for generic rates (where non- 
trivial spatial correlations are expected) are not known. 

In this article we introduce a totally asymmetric ex- 
clusion process on a ring with non conserved internal de- 
grees of freedom. In this model, apart from the usual 
forward hopping, each particle in one of the v possible 
internal states can change to any other v — 1 states with 
different rates. The steady state of the model could be 
calculated exactly. Unlike TASEP on a ring, which has 
an uncorrelated steady state, this model show non-zero 
spatial correlations when hop rate of particles depend on 
their internal states. 



II. THE MODEL 

The model is defined on a one dimensional lattice with 
periodic boundary conditions. The sites, labelled by i = 
1,2 ... L, can either be empty or occupied by at most one 
particle. Each particle can be in one of the v possible 
internal states I = l,2,....v . Correspondingly, at the 
i th site we take the site variable Sj = for vacancy and 
Si = 1,2 . . .v for the different internal states. 

As in TASEP, the particles hop towards the right 
neighbouring site if it is vacant. However, the hop rate 
depends on the internal state of the particle. Correspond- 
ing dynamics is 
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(1) 



where 7 = 1, 2, ....v represents the different internal states 
of the particles and ai are the hop rates of the particles in 
the I th state. In addition to this hopping dynamics, the 
particles are also allowed to change their internal states 
when their right neighbouring sites are occupied. In other 
words 
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J K. 



(2) 



where the rates pu and pji do not depend on the state 
K of the neighbouring particle. Note, that K need not 
necessarily differ from 7 or J. This dynamics @ natu- 
rally do not show up for v = 1. The only dynamics ([1]) 
is identical to that of TASEP. 

For simplicity we first study the case v = 2 in details; 
generalization to arbitrary v is straightforward. Particles 
in the two possible internal states here are denoted by + 
and — for notational convenience; correspondingly the 
site variable Sj = 0, ±. Both + and — particles can move 
only to the rightward vacant site with rates a±. Other- 
wise, if the rightward neighbour is occupied, the particles 
can flip their state (± — > =p) with rates p±. Explicitly, 
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Clearly, the total number of particles N + + N- and 
vacancies Nq are conserved by this dynamics, whereas 
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FIG. 1: Mapping of the exclusion process to the ZRP. + 
and — particles are regarded as two different kinds of boxes: 
shown as square(blue) and semi-circular (red) boxes respec- 
tively. Uninterrupted sequence of 0s in front of ± are consid- 
ered as the particles in the box. Particles in this ZRP can 
hop from the ± boxes to the left only, with rates a±. The 
empty boxes can change state, ± — >• =p with rates p±. 



individually N + and N- are not. Thus, both the den- 



sities fluctuate and only their averages p± 
reach a stationary value in the steady state. 

The model can be mapped to the Zero Range Process 
(ZRP) [Hill by identifying + and — as two different kinds 
of boxes denoted by r = ± respectively. Uninterrupted 
sequence of 0s, say n in number, to the right of a +(— ) 
particle is regarded as a +(— ) box containing n particles. 
The mapping is described schematically in Fig. 1. Thus, 
this ZRP with No particles distributed in M = L — No 
boxes obey the following dynamics: 
(a) particles from a +(— ) box moves to the left box with 
constant rate a+(a_), (b) an empty +(— ) box can alter 
its state with rate p+ (p_ ) . 

Note that here the hop rates a± of ZRP are indepen- 
dent of the number of particles in the departure box. The 
box dynamics (b) here can produce correlations among 
boxes (absent in ordinary ZRP) as the non-empty ± 
boxes can not alter their sign. 



III. STEADY STATE FOR v = 2 

A generic configuration can be written as 
{niTi, 712T2 . . . njwTiw} where n k is the number of 
particles in k th box of type Tk- ft can be shown that 
steady state of this ZRP has a product measure : 

P(n 1 Ti 1 n 2 T 2 , . ..n M T M ) ~ /n ("i)/t 2 ("-2) • • ■ /t m ("m) 

(4) 

where weight of an individual ± box containing n parti- 
cles is f±(n). 



Proof of Product measure 



{. . . nk-iTk-i,rikTk,nk+iTk+i ■ ■ ■}■ Let us construct C" 
for all possible C's. 

Case-I rife = : In this case the k th box can change 
sign. ThusC" = {. . . rc fc _iTfe_i, 0f fc , Uk+iTk+i ...}, where 
Tk = —Tk- A choice C" = C along with Eq. (gj) and © 
gives 



/-(0) 
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P+ 

P- 



(6) 



Case II nfe > : In this case a particle from 
box k can move to k — 1 with rate a Tk . Thus C = 
{... (rik-i + l)r fe _i, (n k - l)r fe , n k+1 r k+1 ...}. We choose 
C" = {. ..n k -iTk-i, (n k -l)T k , (n fc+ i + l)r fc+ i . . . }. Then 
the condition of pairwise balance along with (j4]) demands 



fr k ( n k) 
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fr k+1 (n k+l ) 



(7) 



The constant c is independent of n and r. We set c = a + 
without any loss of generality. Since the steady state 
weights are yet to be normalized, we can take /+(0) = 1 
without loss of generality. Then, Eqs. © and ([7]) result 



/+(n)=l; f-(n)=pa r ' 



(8) 



where a = — and p = — are the ratios of hop rates 
and flip rates respectively. 

Now, weights of every configuration in the lattice can 
be written in terms of the steady state weights of ZRP. 



P({ Si }\i = l,...L) - P({n k Tk}\k = l,...M) 



n k ) 



(9) 



k=l 



Spatial indices of the model Q on a lattice are not 
carried over to their corresponding ZRP version, which 
makes calculation of the spatial correlations unreason- 
ably difficult. We follow a method discussed in [t3 |. 
First, let us rewrite the steady state weights as a ma- 
trix product state, by replacing Si by the corresponding 
matrices A%. . 



P({ Si })=Tr[l[X Si } 



(10) 



Since we already know P({si}) from Eq. ^ in terms of 
f T (n)s, X Si s must be chosen such that 



The proof can be constructed using pairwise balance 
where we find a unique configuration C" for every tran- 
sition C — )• C such that 



P(C)W{C C) = P(C")W(C" -> C). 



(5) 



It is sufficient to consider all possible transitions that 
changes the state n k r k of the box k. Thus C — 
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rr[n^]=n^( n *)- 



k=l 



(11) 



We proceed with a simpler notation X + = D, X- = E 
and Xq = A and the following choice, 



D=\d 1 )(d 2 \; E=\e x )(e 2 \. 



(12) 
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Now, Eq. (fTTj) and (fT2]) together impose the following 
conditions on matrix A, 

(d 2 \A n \d 1 ) = {d 2 \A n \e 1 )=f+(n) = l 
(e 2 \A n \ ei ) = {e 2 \A n \d x ) = /_(n) = P a n (13) 

Note, that unlike the usual Matrix Product Ansatz[(| 0], 
these conditions on D, E and A do not depend on the 
dynamics of the model explicitly. This will be discussed 
in some length later in lTnBl 

An infinite dimensional representation of these matri- 
ces that satisfy Eq. (fT3"j) . similar to those obtained for 
the usual ZRP with one kind of boxes [l2j], can always 
be constructed. Fortunately for this model, as f T (n) are 
simple functions, we have a two dimensional representa- 
tion, 



|di> = |ei) = |l) + |2>;<d a | 
A 



(l|;<e 2 |=p<2|; 
|l)(l| + a|2)(2|. 



(14) 



Here |1) = (J) and |2) = (°) are the standard basis 
vectors in two dimensional vector space. Explicitly, 
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(15) 



B. Matrix Product Ansatz 

In recent years steady state weights of several exclu- 
sion models have been obtained applying Matrix Product 
Ansatz (MPA)@,I3- To use MPA for exclusion models, 
first one writes the steady state weight as the trace of 
the product of matrices, similar to Eq. (fTU|). The ma- 
trices are required to satisfy certain algebraic relations, 
which depends on the dynamics and also involve a set of 
auxiliary matrices. 

In this model with v = 2, we already have matrices 
D, E and A [given by Eq. (fT"5|) ] which provide the exact 
steady state weight of the model in matrix product form. 
However, it is not obvious that these matrices, which are 
constrained by Eq. (ITBl . also follow the algebraic rela- 
tions required by MPA. In the following we show that, 
in fact for this particular dynamics, matrices of (1151) can 
be made to satisfy the required equations with a suitable 
choice of the auxiliary matrices. An explicit represen- 
tation of these auxiliary matrices D, E and A are also 
given. 

The complete set of equations for the matrices D, E 
and A, as required by the MPA, are 



a+DA = 

u-EA = 

p+DD -p_ED = 

p+DE - p-EE = 



-AD + AD = DA- DA 
-AE + AE = EA- EA 
-ED + ED = DD - DD 
-EE + EE = DE - DE, (16) 



unique solution. For simplicity we choose D = D, E = E, 
which reduces the above set of equations to, 



p+DD =p_ED 
(1 - a + )DA = DA 
(l-a-)EA = EA 



p+DE = p-EE 
a+DA — AD = -AD 
a-EA — AE = -AE (17) 



Now we need to check, if D, E and A obtained earlier in 
Eq. (jl5p satisfy above equation along with some auxiliary 
matrix A. 

\l-a+) 

a(l-a!_), 
tently solves the above equation along with Eq. (IT51) . 



It turns out that A 



IV. CORRELATION FUNCTIONS (v = 2) 

To calculate the partition function one needs to take 
care of the conservation of total number of Os. Instead, 
we work in grand canonical ensemble (GCE), where the 
fugacity z associated with 'A's fixes the average density 
of Os. Thus the partition function in GCE is 



Z = Tr[(D + E + zA) L ] = Tr[T L ] 



(18) 



where we have used T = D + E + zA 



1 + z p 
1 p + za 



for convenience. The eigenvalues of T are 

A± = -\l+p+z+az±^(l + p + z + az) 2 — Az{a + p + az)] 

(19) 

rpi r~7 \ T, , \ T. m l l i i l l i - r 




FIG. 2: The main figure compares Eqs. ([21]) and (j24j) 
(solid lines) with {+) (circle) and (H — h) (cross) obtained from 
Monte-carlo simulation of a system of size L — 1000, a — 0.3 
and p = 0.6. The correlation function C++ is shown in the 
inset, where solid line corresponds to Eq. (I25[) . 



where the auxiliaries D,E and A also need to be deter- 
mined. These set of equations do not necessarily have a 
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FIG. 3: The average velocity of 0s, obtained from Eq. (|23[) 
is compared with the numerical simulation for L — 1000, a = 
0.3 and p = 0.6. The inset shows correlation function Coo 
versus po- 
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Calculation of expectation values of the observables are 
straightforward in this formulation. For example, 



Pa 



P+ 



(N ) Tr[zAT 



L-n 



P- 



Z Z 
z[(X+ — aA_) — (1 — a)(p + az)] 

A+(A+-A_) 

Z 

(A+ — az) 



(20) 
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A+(A+-A_) 
The average velocity of 0s is 

v = a-[a(+0) + (-0)] = a- 



(21) 



(22) 



zTr[(aD + E)AT L ~ 2 ] 



Z 

za(X+ + p — az — p(A_ — p — az)) 



Xl(X, 
Two point correlations: 
Tr[DDT L - 2 ] 



(++) = 

C++ = 
(00) = 



A _) 



T L-2 



(23) 



Z 

(A + — az) 
A^(A+-A_) 

(++) - (+) 2 = 

^L-2 



pz(l — a) 
A 2 + (A + -A_)2 



Tr[AAT L - 2 } _ z 2 (T^f 2 + a 2 T^' z ) 



Z Z 
z 2 [(X+ - a 2 A_) - (1 - a 2 )(p + az)} 

A 2 ,(A+-A_) 



(24) 
(25) 

) 

(26) 



Coo = (00)-(0) : 



pz 2 (l — a) 2 

x 2 + (x + -x.y 



(27) 



We have compared these results with those obtained 
from the Monte-carlo simulations. In Fig. 2 we have 
shown variation of (+} and (H — h) with po. The param- 
eters a = 0.3, p = 0.6 and the system size L = 1000 
were fixed. Solid lines correspond to Eqs. (|2"Tj) and (|2"4")l 
respectively, where po was obtained from Eq. (|20p . Inset 
of Fig. 2 compares Eq. (1231) with C ++ obtained from 
simulations. 

Similarly, in Fig. 3 we have compared the average 
velocity of 0s with Eq. (|23p . The inset therein shows 
variation of C 00 = (00) - (0) 2 . 



GENERAL RESULTS FOR v > 2 



In this section we will extend the results of section lTlI Al 
and obtain the exact steady state weights of the model 
with generic v > 2. 

Here, it is evident that the dynamics ([TJ and ([2]) can 
also be mapped to that of a ZRP with v kinds of boxes, 
where a particle from a randomly chosen box of kind / 
moves to the left box with a constant rate ai and an 
empty box can change its internal state. As before we 
look for a factorized steady state satisfying pairwise bal- 
ance condition. Accordingly, for non-empty boxes, we 
obtain 



aj- 



/ /(» + !) 



= c,VJ= 1,2.. 



(28) 



similar to Eq. ([?])• Here again c is an arbitrary constant, 
independent of the number of particles n and the state 
of the box /, can be taken to be 1 without any loss of 
generality. Thus fi(n) — //(0)/a™. The relative weights 
of the empty boxes, again, satisfy 



.MO) 
fj(P) 



P.7I 
PI J 



(29) 



for all possible pairs (/, J) . This condition, which is a 
generalized version of the Eq. ([6]), demands that the 
steady state cannot have a product measure form if all 
the v(y — 1) rates pijs are independent; they must be 
related in the following way 



Pu Pjk _ Pik 

PJI PK.J PKI 



V/, J, K 



(30) 



This set of equations imposes (y — l)(v — 2)/2 con- 
straints leaving (v — \)(yj1 + 1) independent pijs. Fur- 
ther, we take /i (0) = 1 as the weights are not normalized 
yet, which results in 



fi{n) = — and 



fi(n) = — \ V/>1. (31) 

Pn a 1 } 



To hnd the spatial correlation functions we rewrite the 
steady state in a matrix product form which is a direct 
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generalization of the two species case; each particle of 
species I is replaced by D[ and vacancies by A. The 
non-commuting set of matrices Di and A must satify 
Eq. (11) with the above //. This can be achieved by 
choosing D] = \d)(di\, similar to Eq. (|12l) . resulting in 

{d!\A n \d) =/j(n) 



0.02 



which is a generalized form of Eq. (113)) . 

It is straightforward to find a v dimensional represen- 
tation of these matrices; \d) — J^i K)> {di\ = fjf^l an d 
A = Diag( — , . . . —). Here, {\i)} are the standard basis 
set for the //-dimensional vector space. 

Let us discuss the v = 3 case in some details. Explic- 
itly, the dynamics is 



10-^-01; 
1/^2/; 

P21 



20- 
2/4 



30^>03 

P13 



> Q2 02; 

P31 



The model is mapped to a ZRP with three different kinds 
of boxes. Its steady state has product measure only when 
the rates pu follow Eq. (f3"U|) ; 



Pl2 P23 
P21 P32 



Pl3 
P31 ' 



Thus, any five out of these six pijs can be chosen indepe- 
dently and the sixth one is fixed by the above equation. 

The spatial correlations of this system on the lattice 
can be calculated in a straightforward manner using the 
prescription decribed for v = 2. All one needs is the 
explicit representation of the matrices, which are, 



Do 



D, 




A 




(32) 



The correlation functions here are found to be quali- 
tatively same as that of v — 2 case. 



VI. DISCUSSION AND CONCLUSION 

Several multiple species exclusion processes with or 
without conservation have been studied earlier, mostly 
numerically, in different contexts. In this article we study 
a multispecies exclusion model where exact steady state 
weights and spacial correlations are calculated analyti- 
cally. 

This model with v = 1 is identical to TASEP, which, 
on a ring, has spatially uncorrelated steady state. When 
internal degrees of freedom are introduced {%. e. v > 1), 
two new features appear in the dynamics : i) particles 
in different internal states hop with different rates, ii) 
the particles can change their internal states resulting in 




FIG. 4: Spatial correlations of the conserved model obtained 
from Monte-Carlo simulations (symbols) on a system of size 
L — 1000 are compared with exact correlation functions of 
the non-conserved model (line) with p = 0.8 ,13]. For both 
the models a = 0.5. 



non-conservation of number of particles belonging to each 
state. To investigate wheih of these features generates 
spatial correlations we revisit v = 2 case of the model 
with the nonconserving part of the dynamics ([2]) replaced 
by a conserving one, 



±0^>0± 



(33) 



Symmetic exchange of + and — particles ensures that 
particle current gets contribution only from the hoping 
dynamics, as was the case in the non-conserved syetm. 

The most generic particle conserving two species model 
tt' — > t't with rate w TT ', where r ^ t' e 0,1,2 
has been discussed in Section 7 of the review article 
7]. It has been shown there that the steady state of 
this generic model can be obtained using Matrix Prod- 
uct Ansatz, when the rates satisfy a set of conditions, 
Eqs.(7.8)- (7.13) therein. The particle conserving model 
(f3"3"| discussed here is a special case of the above with 
wio = a+,W2o = a_,wi2 = 1 = W21, and all other 
rates are zero. Here, MPA provides an exact solution 
only when a + = a_; corresponding steady state turns 
out to be spatially uncorrelated. This indicates that 
non- conservation is irrelevant, at least for a+ = a_, as 
the non-conserving model ^ too provides an uncorre- 
lated steady state in this case. To understand, if non- 
conservation plays any role in generating spatial correla- 
tions, we study model (|33[) numerically for generic rates 
a±. These Monte-Carlo simulations reveal that non- 
trivial spatial correlations appear for a+ =/= a-. In Fig. [4] 
we compare some of the two-point correlations obtained 
from Monte-Carlo simulations of the conserved system 
with the known exact results of the corresponding [l3| 
non-conserved model. Clearly the spatial correlations of 
the conserved model are qualitatively similar to those of 
the non-conserved case. Thus, conservation of internal 
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degrees is an irrelevant criterion in developing spatial cor- 
relations in exclusion processes, rather the unequal hop 
rates are responsible. 

In conclusion, we have studied an exclusion process on 
a one dimensional system with periodic boundary, where 
particles carry v internal degrees of freedom. Along with 
a directional hopping dynamics that depends on the in- 
ternal state of the particle non-conservation is introduced 
by allowing a particle to change its internal state when 
the target site is occupied. We show that this model 
can be mapped to a zero range process where particles 
distributed in v diffrcnt kinds of boxes hop with different 
rates and empty boxes can change their state. The steady 
state weights of the exclusion process could be written in 



matrix product form using the exact single box weights 
obtained from ZRP. Exact spatial correlations have been 
calculated exploiting the matrix product form. From the 
comparison of these correlations with those of a corre- 
sponding conserved model we conclude that the unequal 
hop rates of particles belonging to different internal states 
is the possible cause for the non-trivial spatial correla- 
tions observed on a ring. 
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